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Abstract 



We extend the effective field theory approach which successfully describes or- 
dinary nuclei and nuclear matter to incorporate strangeness in nuclear struc- 
ture. Central object is a chiral effective Lagrangian involving the baryon 
octet, the Goldstone boson octet, the vector meson octet and a light scalar 
singlet. According to the rules of effective field theory, we include all inter- 
action terms (up to a given order of truncation) that are consistent with the 
underlying symmetries of QCD. We develop a mean-field approximation and 
study nuclear matter as a simple model for multi-strange systems. A D-type 
Yukawa coupling between baryons and vector mesons leads to A — S*^ flavor 
mixing in the nuclear medium. We study flavor oscillations in the nuclear 
matter ground state which are closely related to the phenomenon of neutrino 
oscillations. 
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I. INTRODUCTION 



Strangeness adds another, still largely unexplored, dimension to nuclear structure. On 
the experimental side physics of hypernuclei is approaching a phase in which not only ground 
state energies but also excitation spectra and electromagnetic properties are being mea- 
sured m. To explain properties of hypernuclei, detailed information on the elementary 
nucleon-hyperon and hyperon-hyperon interaction is needed which, at present, is scarce and 
incomplete. 

On the theoretical side strangeness in nuclear structure is studied from a variety of 
different perspectives. Microscopic meson-exchange models have been constructed which 
accurately reproduce the rich nucleon-nucleon and the more scarce hyperon-nucleon data 
IQJ^. Binding energies and single particle spectra of single hypernuclei are described in 
nonrelativistic and relativistic mean-field models. The relativistic models usually 
involve the baryon octet and several strange and nonstrange mesons. A very compelling 
feature of the relativistic approach has been the reproduction of the observed small spin-orbit 
splittings of A hypernuclei by introducing an appropriate hyperon-meson tensor coupling 
M. Another goal of the study of strangeness in nuclear structure is to extrapolate to multi- 



strange systems |]10[- The existence of a very large class of bound, multi-strange objects has 



been suggested which might be created in central collisions of very heavy ions |]TT 



At present, double A hypernuclei are the only source of information on multi-strange 
systems. Only a few events have been identified indicating a strong attractive AA inter- 
action. While mean field models well reproduce properties of single A hypernuclei, theoretical 
uncertainties arise in extrapolations from single to multi-strange systems. The respective 
information on the AA interaction derived from the few double A hypernuclei appears to be 
rather ambiguous. For example, nonrelativistic potential models |jT2| reproduce the binding 
energies with AA potentials of completely different types. 

Neutron stars are another area where the study of strangeness has received considerable 
attention due to the possibility of kaon condensation |T^. Uncertainties arise from an 
imperfect knowledge of the equation of state when hyperons are present [14|. Particularly, 
the question of kaon condensation is very sensitive to specific model features [T^JT^ related 
to the hyperon-nucleon and hyperon-hyperon interaction. 

To summarize, a more systematic theoretical treatment of strange systems is needed to 
achieve more predictive power and to guide present and future experiments. 

The concepts and methods of effective field theory (EFT) successfully describe the low- 
energy phenomenology of quantum chromodynamics (QCD). Motivated by this success, EFT 
concepts have recently been applied to models of nuclear structure [jl^. A chiral effective 
Lagrangian for ordinary nuclei has lead to new insights into models based on Quantum 
Hadrodynamics (QHD) |jl8|. One important implication is that all interaction terms that 
are consistent with the underlying symmetries of QCD should be included. 

An EFT for strange systems is more involved due to the vast number of allowed inter- 
action terms arising form the underlying SU{3) group structure. However, it is unnatural 



for some interaction terms to vanish without a relevant symmetry argument |]T9|. From 
this modern point of view most relativistic mean-field models which have been employed 
to describe strange nuclear systems (for instance P,|6|, p!0|j20[] ) are incomplete because only a 
very restricted subset of the allowed interaction terms are considered. 
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It is our goal to extend the ETF formulated in Ref. JT^ to describe systems with 



strangeness. The hadronic degrees of freedom are taken to be the baryon octet, the Goldstone 
boson octet and the vector meson nonet. In addition a light scalar field is included which 
simulates the exchange of correlated pions and kaons. In an EFT approach the structure of 
the particles is described with increasing detail by including more and more interactions in 
a derivative expansion. This implies that the EFT will generally contain an infinite number 
of interaction terms, and one needs an organizing principle to make sensible predictions. 
Applications to ordinary nuclei ||21[] have demonstrated that the meson mean fields provide 
useful expansion parameters which allows a truncation of the EFT. 

A framework which includes the most general types of interactions leads to new and 
interesting many-body effects. Most prominently, the D-type coupling between baryons and 
vector mesons gives rise to A - S° flavor mixing. Although it arises naturally in the EFT 
description flavor mixing has never been studied in this context. 

To provide a first orientation of the EFT approach we will study strange nuclear matter 
as a simple model for multi-strange systems. The central point in the discussion is the 
analysis of A — flavor mixing. The physical nature of this effect is similar to the recently 
much discussed neutrino oscillations ||2^. The primary result is that nuclear matter is 
generally in a state of mixed flavor rather than in a state with distinct A and T,^ particles. 
A particle interpretation is only possible in terms of the actual mass eigenstates which are a 
superposition of the flavor eigenstates. As a consequence, systems which contain A hyperons 
always have a small admixture of S° hyperons. Moreover, disturbing the time independent 
nuclear matter ground state leads to flavor oscillations characterized by distinct frequencies. 
Flavor mixing is driven by the mean field of the p meson and does not occur in isospin 
saturated systems. At low density the effect is very small, but we expect that signatures 
of this new feature will survive in heavy and very asymmetric hypernuclei. For instance, 
A — S° flavor mixing could produce considerable deviations of the magnetic moments from 
the Schmidt values 



The idea of A — S mixing is not new. In the vacuum it arises from small isospin viola- 
tion induced by the mass difference between up and down quarks |24] and electromagnetic 
interactions I^S] which account for electromagnetic mass differences between members of 
the same isospin multiplet. Hence, explicit 5'f/(3)-symmetry breaking is responsible for the 
mixing. In contrast, the effect we study is a true many-body effect which also arises in the 
chiral limit. Moreover, the effect is considerably larger than the vacuum mixing. 

A systematic treatment of strangeness based on an EFT introduces a large number of 
coupling constants. For the parameters in the nucleon sector we rely on the studies in 
Ref. |]T^. Ultimately, the coupling constants in the hyperon sector have to be constrained 
by properties of hypernuclei. Work along these lines has been reported in Ref. ||20|. The 



couplings of the A hyperon were adjusted by a least square fit to experimental single particle 
spectra of selected hypernuclei. However, these couplings cannot be used in our framework 
because the analysis was based on a model which includes a very restricted subsets of the 
allowed interaction terms. We therefore resort to a more phenomenological approach which 
is based on the hyperon potential in nuclear matter. Experience has shown that different 
models which give a realistic description of single hypernuclei predict very similar values 



for the hyperon potential in nuclear matter |26,27|. Although this simple phenomenological 



approach cannot constrain all the relevant parameters, particularly the nonlinear meson- 
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meson couplings, the nuclear matter studies provide a useful qualitative description of the 
pertinent many-body effects. 

More recently, an approach similar to ours has been proposed based on linear [|^ and non- 
linear P] realization of chiral symmetry combined with the idea of broken scale invariance. 
Although a framework which incorporates broken scale invariance successfully describes 



properties of finite nuclei pl| , |28[| the concept appears to be less compelling p9[. Compared 
to chiral symmetry, breaking of scale invariance is much larger on the scales relevant in 
nuclear structure. Symmetry breaking terms are essential for predicting masses and cou- 
pling constants and it is unclear to which extent the symmetry pattern survive at the end. 
Moreover, we will argue that even for the case of approximate chiral symmetry important 
parts of the effective theory, namely the nonlinear meson-meson interactions, appear to be 
rather asymmetric. The relevant couplings receive contributions from a very large number 
of terms in the original Lagrangian. If these couplings can be constrained by reproducing 
properties of normal and hypernuclei, complicated many-body effects will be included in an 
approximate form and it will be very difficult and also unnecessary to disentangle how these 
couplings arise form the underlying symmetries. 

The outline of this paper is as follows: In Sec. 0, we present the effective Lagrangian 
which lies at the heart of the EFT. Section ^ contains the mean field approximation. In 
part Llllj we focus on the derivation of the thermodynamic potential in the presence of 
A — S'^ flavor mixing. In part |III Q we present the relations that determine the equation of 
state. In Sec. |^ we discuss how the model parameters are specified. Section deals with 
strange nuclear matter as a simple model for multi-strange systems; various systems with 
different isospin and strangeness content are discussed, and the impact of the flavor mixing 
on the flavor content is illustrated. In Sec. 0, we apply our model to describe dense matter 
in neutron stars. In Sec. |V11| we study flavor oscillations in nuclear matter. Section |V111 
contains a short summary. 



II. THE EFFECTIVE LAGRANGIAN 

Recently, it has been shown that hadronic phenomenology can be combined with the 
ideas of EFT in form of an effective Lagrangian that realizes chiral symmetry and vector 
meson dominance p!7|JT8[] . This approach has been successful in describing the properties 



of ordinary nuclei and nuclear matter. In the following we will generalize these ideas to 
describe strangeness in nuclear systems. Although some of the material with restriction to 



the nucleon sector can be found in Refs. |17,18|, we repeat it here to develop new aspects 



arising from the more involved SU (3) structure. 

The effective degrees of freedom in our approach are the baryon octet, the Goldstone 
boson octet and the vector meson nonet. In addition we also include a light scalar field 
to simulate the exchange of correlated pions and kaons. In principle, the dynamics of 
the non-Goldstone bosons could be generated through pion and kaon loops. In studying 
the many-body problem it is of advantage if the evaluation of complicated loop integrals 
can be avoided. Moreover, the midrange part of the baryon-baryon interaction is efficiently 
described by the exchange of vector mesons and a light scalar meson. In the meson sector this 
picture is confirmed by the observation that the coupling constants can be understood from 
vector meson exchange pO] . Constraints are imposed through the underlying symmetries of 
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QCD including Lorentz invariance, parity conservation, and approximate chiral symmetry; 
and all allowed interaction terms must be included. To write down a general effective theory 
an organizing and truncation scheme is needed for the theory to have any predictive power. 



Following Ref. we assign to each term in the effective Lagrangian an index 



d+- + b 



(1) 



where d denotes the number of derivatives, n is the number of baryon fields, and b the number 
of non-Goldstone boson fields. Derivatives acting on baryon fields are not included in d 
because they generate powers of the large baryon masses, which spoil the counting scheme. 
In principle, an infinite number of terms is possible and a meaningful way to truncate the 
effective theory is needed. It has been shown that at low and moderate nuclear densities 
the meson mean fields and their gradients are sufficiently small to provide useful expansion 
parameters If one assumes that the coefficients of any term in the Lagrangian are 

natural, the theory can be organized in powers of the fields and their derivatives. Truncation 
of the Lagrangian at a given order in u leads to a finite number of parameters which have 
to be determined by nuclear observables. Furthermore, if naturalness holds, the omitted 
higher order terms are small. 

The effective degrees of freedom are introduced as nonliner realizations of the chiral 
group SU (3) L X SU{3)ji. For the Goldstone bosons we introduce 



u 



(2) 



where 11 is a traceless hermitian 3x3 matrix in flavor space 



n 



/ J_„ + J_7rO 

7r~ 



vr 



V 



K- 



(3) 



The baryons are collected in a 3 x 3 traceless hermitian matrix B 



B 



V 



1 vo 
73^ 



P 

n 



7-.^ -72^" 
2° 



(4) 



The transformation properties of the Goldstone bosons and the baryons under the chiral 
group are |^ 

qlu h\li, g) = h{U, g)ugli, 
B^hiU,g)Bh\U,g) , 



(5) 



with {gi, gji) G SU (3) l x SU{3) r and h{Il, g) is the so-called compensator field representing 
an element of the conserved subgroup S'f/(3)y. For the vector meson octet 



/ J_T/8 I J_^i 



P 



K*+ 
Kf 



J_V^8 _ J_ 

K*- IC*^ "^V^ I 



(6) 
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we assume the standard transformation properties of any matter field 

V,^h{Yi,g)V,h\n,g) . (7) 

We also introduce a vector meson singlet and a light isoscalar scalar meson (p. The 
physical uo and mesons arise from the mixing relation 

uj^ = comS^ + ^m{e)V^ , (8) 
0^ = sin(^)^^ - cos(^)V;^ . 

The effective Lagrangian can be decomposed in two parts: 

C = Cb + Cm , (9) 

each part will contain terms up to a maximum value z/ = 4. The part involving the baryons 
is given by 



Cb = Tr 



B{zl^ - Mo)b] + «pTr(s^75[4, B]) + aoTr^Bt-f^i^, B}) (10) 
gFTT(B[Y,B]) - goTT(B{Y,B}) - gsTiiB^B) 



+ 5C 



SB 



Most terms in this part of the Lagrangian are standard |3^. The covariant derivative is 
defined by 

D^,B = d,,B + [r^,B] , (11) 

with the connection 

T^ = ^iu^d,,u + ud^u^) . (12) 

In addition to the covariant derivative the baryons couple to the Goldstone bosons via 

= ^{u^df,u - ud^v)) . (13) 

The symmetry breaking part of Eq. (ID) contains the quark mass matrix 

M = diag{m„,mrf,mj (14) 

and generates the baryon masses and scalar couplings 

Tr(BMoB) - dCl"" = N{Mm - g%v)N + A(Ma - (7»A (15) 

+ S(M2 - g^ip) S + H(Ms - g^cp)E . 

The couplings to the scalar field are purely phenomenological simulating the rather complex 
process of correlated two pion and two kaon exchange. The couplings gp are not constrained 
by chiral symmetry and have to be determined for each baryon flavor separately. All the 
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terms in Eq. ([T5|) can be constructed by using suitable combinations of the baryon matrix 
and the quark mass matrix. For the scalar couplings terms quadratic in Ai are needed, i. e. 



9f = 0{ml) . (16) 

The terms listed in Eq. ([T0| ) are all of order u < 2. The ellipsis stands for terms which 
contain one or more derivatives of the meson fields, e.g. tensor couplings, and for the 



couplings to the electromagnetic field. Although important for finite nuclei calculations fl? 
these terms are not needed for our discussion of nuclear matter. Also omitted are terms 
which contain higher-order couplings between Goldstone bosons and contact interaction of 
four and more baryons. Contact interactions are taken into account by the couplings to the 
scalar and vector mesons |r7[ . 

The mesonic part of the Lagrangian can be written as 

Cm = -If^Tt{{A„ A^) - ^Tr( V^'^^ - 2mlV,Vn - \s,.S^'' + ImlS.S^ (17) 
1 
2 



+ -{d,^d^^ - my) + sell + V(¥', V„ S,) 



with 



V^. = D^V,-D,V^ and S^, = d^S, - d,S^ , (18) 



for the abelian field tensors. Explicit symmetry breaking terms are collected in SCff^ which, 
in analogy to the baryons, generate the physical meson masses. For the vector mesons this 
leads to 

^mlTT{V,Vn + ImlS.S^ + 6Cff^ (19) 

Nonlinear meson-meson interactions are collected in the potential 

V(^, V,, S,) = -1/^3^' - ^K4/ (20) 

+ ^US.Sn' + |c8(Tr(V,V'^))^ + \,,TT{y,V^)S,S^ 

+ l=^,i:i{v,v^v,)s^ + . . . + . 

The part 5Cff^ contains additional symmetry breaking terms with three and four meson 
fields and will be discussed in Section 0. The list of terms in Eq. ( pOD is not complete. The 
ellipsis stands for additional terms with four octet fields and terms which involve gradients of 
the meson fields. The gradient terms can be disregarded in the discussion of nuclear matter. 
Furthermore, only the time-like component of the meson mean fields is nonvanishing and 
the omitted contributions involving four octet fields can be reduced to the given terms. 
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III. THE MEAN FIELD APPROXIMATION 



A. The loop expansion 

In the following we will disregard weak decays and small isospin violations. We con- 
sider nuclear matter as a system with conserved baryon number (Nb), isospin (T3) and 
strangeness [S). The corresponding chemical potentials are introduced by adding to the 
effective Lagrangian in Eq. (^) the contribution 

HbTt{WB) + ^3Tr(S7°[T3, B]) + fisTT{B^\B - [Y, B])) , (21) 

where S = Nb — Y was used. The isospin and hypercharge operator are expressed in terms 
of Gell-Mann matrices 

n = lx3 , Y = ^Xs. (22) 

The mean-field approximation follows from the the one-loop contribution to the thermody- 
namic potential. For the main parts the derivation is straightforward and we will only 



briefly sketch the steps. Starting point for the loop expansion is the effective Lagrangian in 
Eq. (^ with all the meson field operators shifted by their expectation values: 

n^fi + n , v^^v^ + v,, 

Sf,-^ Sf, + and + (f , 

where the field operators are denoted with a hat. The symmetries of infinite nuclear matter 
simplify the discussion considerably. Translation and rotational invariance demand that the 
expectation values, or mean fields, of all three-vector fields vanish. The mean fields of the 
Goldstone bosons vanish because the nuclear matter ground state is assumed to have good 
parity. Finally, since the third component of the isospin and the total strangeness is fixed 
via the corresponding chemical potentials, SU{3) symmetry demands that only the time- 
like components of the neutral vector mesons have a nonvanishing expectation value. At the 
one-loop level the thermodynamic potential is then obtained by diagonalizing 

£W = Tr (bI^z^B + fiBl^'B + /i37°[T3, B] + f,sl'{B - [Y, 5]))^ (23) 

- gFTT{W[Vo, B]) - goTT[W{Vo, B}) - gsTT{WSoB) 

- N{Mn - g%^)N - A(Ma - glv)^ - S(Ms - ^?|,<^)S - H(Mh - ^4¥^)S 

+ if^l^l + if^l^fl + if^lpl - if^W + V(</?, UJQ, 00, Po) , 



with 



Vo 



^^2 + % 

10 PO. Q 

n 2_T/8 



(24) 
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The corresponding mean fields for tlie uj and arise from the mixing relation Eq. (P) 



CUo = cos(0)5o + s\Yi{e)V^ 
00 = sin(^)S'o - cos(^)l^o^ 



(25) 



The Lagrangian £'^1 is diagonal in flavor space except for a term which mixes the A and TP . 
For the pure flavors the calculation is straightforward and the thermodynamic potential can 
be written as 



VL _Q.p VLn ^T.+ ^s- ^H- ^ASO 



(26) 



m, 



- 7J"^0</'O - T^T^pPo + T^^sV - V(<^, UJQ, 00, Po) 



The one-body contribution can be divided into an explicit density dependent part and a 
divergent vacuum part 



with 



I J- , l™*^ 

Up + Up 



(27) 



u\ 



Stt^ jm 



up 



dE{E^ 



m;2)3/2 



(28) 



and 



u. 



1 r(2-f) 

27r2 D{D - 2) 



* 4 





D/2-2 


-Ml, 




D/2-2\ 






[47rA2j 





(29) 



The effective chemical potentials and the effective mass for each individual flavor are listed 
in Table |I[ The vacuum contribution includes a vacuum subtraction and was calculated 
in dimensional regularization which introduces the mass parameter A. Before we discuss 
Eq. ( P^D further let us turn to the contribution of the A and S*^ which is more complicated 
due to the flavor mixing. 



B. A — S*^ flavor mixing 

Flavor mixing results from the nondiagonal part of the one-loop Lagrangian Eq. (p^) 
'[1] _ vt7. _ ^,0t/0 _ i\/r*\-a, . i \n^,J^i^ _ ^,0t/0,. 



£W = ^A(^^ - 7V^° - MX)^A + ^so(^^ - YVio - M^)^so (30) 
with 

= gi^p' , (31) 

and where the other potentials and effective masses are listed in Table |. At the one-loop 
level the mean fields are the baryon self energies arising from a resummation of tadpole 
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diagrams involving loops with one baryon propagator. Due to the coupling (^^^ the vector 
self energies are non-diagonal in the A — sector of flavor space. 

The Lagrangian Eq. ( pOD leads to a system of coupled Dirac equations 

- /y^o - Ml)vl/^ = 7V^vl/^o , (32) 
- 7Vs°o - M*)*so = 7V°vl/A , (33) 

which has four energy eigenvalues corresponding to two particle and anti-particle solutions. 
Each flavor is then represented as a superposition of these two solutions. To be more specific 
let us consider the chiral limit M\ = = M*. In this case the algebraic structure simplifies 
considerably and the Lagrangian Eq. (|30|) can be fully diagonalized by substituting for the 
fields 

= cos(a)^i + sin(a)^2 (34) 
= - sin(a)^i + cos(a)^2 (35) 
4Lo = - - M*)^i + ^2(i^ - 7°^2° - M*)^2 . (36) 

The mixing angle is determined by 

tan2(a) - ^ ^° tan(a) -1 = 0, (37) 



' m 



and the potentials are given by 



= cos\a)V^ + sm\a)V^o - 2 sin(a) cos(a)K° , (38) 
= sm'^{a)V^ + cos\a)Vio + 2 sin(a) cos(a)V;^ . (39) 

By using plane wave solutions for the corresponding Dirac equation the energy eigenvalues 
for particle 1 and 2 are now straightforward 



= y° ± Jp2 + jvf*2 , = ± + M*2 . (40) 



The A — S'' mixing is related to the phenomenon of neutrino flavor mixing which gives rise 
to the recently much discussed neutrino oscillations. However, the origin of both effects is 
fundamentally different. Neutrino oscillations are assumed to occur in the vacuum arising 
from a nondiagonal mass matrix in flavor space which contains the vacuum mass parameters 
P2|| . This effect can be appreciably enhanced when neutrinos pass through dense matter as 
predicted by the MSW effect H. In this case the mixing transformation which corresponds 



to Eqs. (pl) and (|35|) induces a nontrivial vacuum structure with two unitary inequivalent 



vacuum states corresponding to the flavor and mass representation ||35|. In contrast, the 
A — S° mixing is a true many body effect arising from a nondiagonal vector self energy which 
is generated in the medium. As long as small isospin violations can be neglected, the vacuum 
self energies are diagonal in flavor space and the asymptotic states can be properly identified 
as the pure flavor states. 



As a consequence of the mixing transformation in Eqs. (133) and (PH) the individual flavor 



states arise as superpositions of the particles 1 and 2. This implies that the nuclear matter 
ground state is generally in a state of mixed flavor rather than in a state with distinct A and 
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S*^ particles. A (quasi) particle interpretation is only possible in terms of the actual mass 
eigenstates 1 and 2. 

For the general case we could not find a unitary transformation which leaves the kinetic 
terms in Eq. (|30|) invariant and which decouples the Lagrangian in a form similar to Eq. (p6|). 
However, it is possible to solve the coupled Dirac equations Eqs. (|3^ ) and ( p3| ) and to 
compute the thermodynamic potential. This is most easily done in euclidean space by 
using a path integral representation. The thermodynamic potential is then given by the 
determinant of the (euclidean) Lagrangian in Eq. (pOl). We find 

^ = -2 In ([(/ + if/°)^ +f + MX'] [(/ - i[/°)^ +f + M* ^] (41) 

- V^''[2MIM^ + 2/ - \/^2 _ 2?7f - 2/2]^ - VEV , 

where we have introduced 

/ = + if/0 , 

and where VEV represents a vacuum subtraction. The momentum integral in Eq. (^Tj) is 
performed in euclidean space. The matter and vacuum contribution can be extracted by 
deforming the integration contour down to the real axis, i.e. by using 



/ 1T-ln[/(po-i^)] = W IT- In -TT— +/ — ln[/(po)], 42) 

where we have assumed that f{ipo) is analytic off the real axis and that > 0. Applied to 
Eq. (|1]) this leads to 

^ = u;l^o+u;T^o. (43) 

The density dependent part can be written as 

d^p 



^!(.o = -2 E l^,{U', - Et)Q{U', - Et) , (44) 



where E^ are the eigenvalues of the particle solutions of the Dirac equations Eq. (|32D and 
Eq. ( |55D . These eigenvalues are the roots of a fourth order polynomial given by the argument 
of the logarithm in Eq. (|iTD. The roots have no simple analytic expression but can be 
calculated numerically. Including the vacuum subtraction the vacuum part is given by 

d'^p 



^AEO = -2 /^|ln(^[(/ + iUlf+i + - -^U'^Y + i + M^^] (45) 

- V^'[2MIM*^ + 2/ - _ 2f;02 _ 2p02j^ 

- In f [(/ +^ + Ma^] \p'+p^ + Ms^]^ I . 
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To study the vacuum contribution in more detail it is useful to separate the contribution 
which does not depend on the vector potentials 



with 



cuXso =c^i+^^2 , (46) 



and 



ln(^[(p° + iU^f+f + - '^U^y+f + M^'] (48) 

- V^2[2MXM* + 2p2 - 1/^2 - 2f/f - 2p°2]^ 
- In f +p' + Ml'] [p° + p' + M* I . 



The contribution ui is identical to the result in Eq. (p9D for the pure flavors. The second term 
Eq. (^) vanishes in the chiral limit = and for Vm = 0. In the vacuum uj2 generates 
the n-point functions of the vector mesons at zero four momentum. These functions are 
in general non-zero because the vector mesons couple to a non-conserved current. As a 
consequence transversality is lost. To extract the divergencies Eg . (^81) can be expanded in 
powers of Vj^. Only the first term in this expansion is divergent 

^T^ = ^r(2 - ^)V^{Ml - M^f . 

indicating a renormalization of the p-meson mass. 

From the explicit form of the thermodynamic potential in Eq. (|^) various densities can 
be calculated by taking partial derivatives. The baryon densities for the A and TP follow 
from 

Pb = -T7^ ' Pb = -77^ • (49) 

The corresponding scalar densities are given by 

^ ^ 1 dVtAT.0 A ^ 1 dVtAT.0 

~ V dMl ' ~ V dM^ ■ ^ ^ 

As a consequence of the flavor mixing the mixed A — S'^ baryon density is nonzero: 

pr G|^A7°^^|G > + < G|^°7°^a|G >= . (51) 

In the limit = = M* the flavor densities can be related to the densities of the mass 



eigenstates by using the relations Eq. ( plD and Eq. 
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Pb = cos^(a)p^ + ?,iv?{a)p\ , (52) 
p| = ?,m^{a)p\ + cos^(a;)p| , (53) 
= -2 sin(a) cos(a) [p], - p|) , (54) 



with 



pb-v dv^ ■ ^^^^ 

Thus, the sum of particles 1 and 2 equals the number of A and flavors in the system 

Pl = pU pf ■ (56) 
It is important to note that this relation also holds in the genral case ^ M^. 

C. The equation of state 

Before we display the relations which determine the equation of state let us briefly discuss 
the role of the vacuum contribution. The divergencies can be eliminated by introducing suit- 
able counterterms in the original Lagrangian Eq. (^). At the one-loop level the counterterms 
are a fourth-order polynomial in the scalar field and the p-meson field. In practice, however, 
an explicit calculation of the counterterms is unnecessary. This is based on the observation 
that the vacuum contribution can be expanded in powers of the meson mean-fields. The 
original Lagrangian contains only terms of the order u = 4 and it is therefore consistent to 
truncate the expansion at fourth order. The truncated vacuum part can then be combined 
with the tree level contributions by redefining the coefficients of the nonlinear potential V 
in Eq. d^), i.e. 

The details of how these coefficients arise in the formalism are unimportant because at the 
end of the calculation the coefficients are determined by fits to nuclear observables. 

The actual mean-field configuration is determined by extremization of Eq. (|26| ) at fixed 
chemical potentials. This leads to the self-consistency equations 



rrippo + 



9%ip'B + Pi) + 91pb + g^ipf + pf 


+ pf 


) + gsipf + pf) , 


(57) 


19M-Pb)+9'ApT -Pf) + l 


gki.pt 


-pD+^asPb''" , 


(58) 


gUp'B + Pi) + gtp'^B + 9Upf + pI~ 


+ pT 


) + gtipf + Pb') , 


(59) 


grip's + p:) + glps + ghipT + pT 


+pr) 


+ gsipf + pT) ■ 


(60) 



-'i^-i 

dpo 

d(f)0 

2 dV 

For the pure fiavors the densities on the right hand side are defined by: 
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for the baryon densities and 

F Ml f'F 



Pl = ^,{-1-Mffl\ (61) 



dEiE^ - M;2)^/2 ^ (g2) 



for the scalar densities. For the A and S° the corresponding quantities are hsted in Sec- 
tion [111 B| in Eqs. (^9D-(|5TD. Note that the p-meson couples to the mixed density which 



enters the right hand side of Eq. (|58|) . The self-consistency equations Eqs. (|57D -(|60|) together 
with the expression for the thermodynamic potential Eq. (|26| ) allows the computation of all 
thermodynamic quantities. For example, the energy density follows from the thermodynamic 
relation 



^ = ^ + PbPb + yUsPs + PsPs , (63) 



with the total baryon density 



the total isospin density 



and the total strangeness density 



Pb = Y.pI^ (64) 

F 



pl = T.t!pB, (65) 



pI = J:^^Pb, (66) 

F 

where was introduced for the isospin and for the strangeness characterizing each flavor. 



IV. MODEL PARAMETERS 

Although a Lagrangian is at the heart of the effective hadronic field theory, an expansion 
in powers of the mean fields is a low density expansion, and neglecting many-body effects and 
loops involving the Goldstone bosons is hard to justify. The success of relativist ic mean- 
field models can be understood in the context of density functional theory ||36|. Central 



object is an energy functional of scalar and vector densities. Extremization of the functional 
gives rise to Dirac equations for the baryons with local scalar and vector potentials, not 
only at the one-loop level, but in the general case as well. In the so-called Kohn-Sham ||37 



approach to density functional theory one introduces auxiliary variables corresponding to 
the local potentials. The exact energy functional has kinetic-energy and Hartree parts which 
correspond to the one-loop contributions derived in section |T|, plus an exchange-correlation 
functional which contains all the complicated many-body effects. Formally, the solution 
of the many-body problem can be cast into solving the noninteracting problem of baryons 
moving in the local Kohn-Sham potentials. The resulting Dirac equations have the same 
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form as in a mean-field calculation, but correlation effects can be included if the proper 
exchange-correlation functional can be found. In the effective hadronic theory the meson 
mean fields play the role of Kohn-Sham potentials and by introducing nonlinear meson- 
meson couplings one can implicitly include additional density dependence. Thus, rather 
than try to construct an energy functional from an underlying Lagrangian, the basic idea is 



to approximate the functional using an expansion in terms of the meson mean fields p8| , |39 
If the parameters can be fit to nuclear observables, complicated many-body effects arising 
from loops will be incorporated. 

Guided by these general remarks we will now specify the parameters contained in the 
basic Lagrangian Eq. (^. The baryon and vector meson masses are taken to have their 
experimental values which are listed in Table |T|. The assumption of SU (3) symmetry implies 
that the couplings of the vector mesons to the baryons are characterized by four parameters: 
the octet couplings Qf^Qdi the singlet coupling gs and the mixing angle 6 which relates the 
physical vector mesons to their pure octet and singlet counterparts. The corresponding 



relations are listed in Table |ITT[ For the mixing angle we used the empirical value 9 ~ 37.2° 
P0| . The parameters gPidD and gs were chosen to reproduce the values for the NNuj and 
NN p couplings which have been obtained in Ref. [^] based on an effective field theory for 
ordinary nuclei. As a third constraint the OZI rule is implemented by requiring that the 
nucleon coupling to the (f) meson vanishes. The couplings of the other baryons then follow 



from the SU (3) relations in Table |T] 



In the scalar sector only the scalar coupling of the nucleon and the mass of the scalar 
meson have been determined within the framework of effective field theory [0 . In principle, 
the scalar coupling of the A has to be constraint to reproduce properties of hypernuclei. Such 



work has been reported in Ref. ||20|, however, their couplings cannot be used in our framework 
because the fits were based on a model which includes a very restricted subset of nonlinear 
meson-meson interactions. We therefore resort to a more phenomenological approach 



which will also lead to an estimate for the corresponding couplings of the E and S hyperons. 
There is a considerable amount of data available on binding energies and single particle 
levels of A hypernuclei ^ which are successfully described in various mean-field models. 
The key observation is that different models predict values for the potential felt by a single 
A in nuclear matter within a fairly narrow range [pGj^ : 

Uf, = glif - g%uj^ ^ 27- 28 (27.5)MeV . (67) 

Although the experimental status on S hypernuclei is still controversial, studies of level 
shifts and widths of S~ atoms suggest |^ 



Uj^ = g'^ip _ g^oj^ ^ 20 - 30 (25)MeV . (68) 

The few events which have been attributed to the formation of a H~-hypernucleus can be 
interpreted in terms of a potential |^2| 



= gt^ip - g!^uj° ^ 20 - 25 (22.5)MeV . (69) 

The mean fields and ip entering the potentials are taken at nuclear matter equilibrium and 
depend only on the nucleon couplings. Thus, the hyperon scalar couplings are completely 
determined by Eqs. (|67D-(|69|). For the actual calibration the values in brackets were used. 
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Let us now turn to the nonlinear meson-meson couplings which are collected in the 
nonlinear potential V = V{ip,ujQ,(j)Q, po). In the effective field theory for normal nuclei 
(00 = 0) this potential was determined in a form 



V{(f, Uq, 00 = 0, Po) = ~^^3^^ - ^'«4¥'^ (70) 



11 11 

/I 2 2\ 2 /- 4- 1 2 

+ + )^o + ^C^o + 2^pVPo 



The other allowed terms 



22 22 1 4 

^ Po > ^ Po and , 



can be expected to give negligible contributions in ordinary nuclei and were not included in 
the analysis. In our calibration procedure we require that the nonlinear potential in Eq. (pO[) 
reduces to the form given above for 0o = 0. This determines nine parameters from the set 

{k3, K4, 7]l, r/2, r]i, vl, Ci, Cs, ?2, ^3} • 

We have arbitrarily chosen to let ^3 undetermined. In the general case, (0o 7^ 0), this 
procedure leads to predictions for the various couplings involving the meson based on 
5*^7(3). Complications arise from terms which are linear in the meson mean field, e.g. 
terms of the form 

00 u;^ , 00 etc. 

These terms lead to abnormal, nonvanishing solutions for the mean field 0o if the total 
strangeness density of the system is zero. This is an indication that symmetry breaking 
terms are important for determining the nonlinear potential. The troublesome terms can be 
eliminated by constructing contributions which contain the quark mass matrix and combi- 
nations of three and four meson fields leading to a large number of new and unconstrained 
parameters. As a minimal approach we simply subtract the terms linear in the meson 
field. This leads to 

4 1111 

V(y^,a;,0,p)= CM>,/77</'V^'7:y0'7^P'' (71) 

with i + j + A; + 2/ = 3, 4 and k,2l>2. 

Ultimately, the coefficients Cij^k,i have to be constraint by properties of normal and hyper 
nuclei. The coefficients Cij^k,i arise from a very large number of couplings in the original 
Lagrangian and are essentially independent. For nuclear structure calculations, however, 
it is not necessary to know how these coefficients arise. It is sufficient to know that the 
potential rewritten in terms of the physical meson fields is of the form given in Eq. ([n]). 

The analysis of finite nuclei in Ref. [|17| lead to two sets, Gl and G2, for the parameters 
{k3, K4, r/j^, i]^, C^, r/^} [| which are listed together with the scalar mass and the nucleon- 



meson couplings in Table IV. The parameters were determined by calculating a set of 



The parameter sets Gl and G2 contain aditional parameters involving tensor couplings and 
meson-field gradients which are not needed in nuclear matter. 
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nuclear properties for a selected set of nuclei and by adjusting the parameters to minimize a 
generalized x^- Among others, the nuclear properties include binding energies, rms charge 
radii and spin-orbit splittings. In spite of the differences in the parameters the sets Gl and 
G2 yield similar properties of nuclear matter. 

These sets Gl and G2 are the input in the calibration procedure as outlined above. 
The symmetry breaking parameter ^3 which is not determined is set to zero. For the octet 
couplings we find gp/igp + Qd) = 0.97 and gp/igp + Qd) = 0.95 for the set Gl and G2 
respectively, which is close to the universality value gp/ {gp + gp,) = 1 To test whether 
the calibration procedure leads to reasonable values for the nonlinear coupling constants we 
have checked naturalness. A generic term of the potential can be written as 

The overall coupling g is dimensionless and of 0{1) if naturalness holds. For normal nuclei 
it has been demonstrated that the parameter set Gl and G2 are natural |]T^. We found 
that the additional parameters which contain the meson field are also of 0{1) and that 
our calibration procedure gives (at least) a natural set of couphng constants. 



V. STRANGE NUCLEAR MATTER 

We consider nuclear matter as a simple model for multi-strange systems which might be 



created during the evaporation process following central collisions of very heavy ions |]TT 
We assume that the time scale during which the system is observed is small enough so that 
the weak decays of the strange baryons can be neglected. The system is then characterized 
by the overall strangeness and isospin. To be specific we introduce the ratios 

X3 = ^ = ^ and = (73) 

Nb Pb Nb Pb 

where the densities were introduced at the end of sect. 4, and assume that the isospin 
and strangeness ratio is held constant. Typical values for the isospin ratio in normal and 
hypernuclei are \x^\ ^ 0.11. Strangeness ratios in heavy hypernuclei are very small but values 
up to 0:5 ~ 1/3 have been observed in light single and double hypernuclei |jl|. Theoretically, 
metastable multi-strange systems have been predicted with strangeness ratios exceeding 
xs = 1 [0. 

At this point some caveats must be added. Our discussion of the theory will encompass 
a wide range of X3 and xs and we will extrapolate to regimes of high densities to obtain 
estimates for the empirical size of the new effects. The model was calibrated by using 
information of normal nuclei and single hypernuclei. Terms which have been neglected or 
which cannot be calibrated accurately, i.e., nonlinear couplings involving the p and meson, 
are likely to be more important in systems with large isospin and strangeness ratios than in 
ordinary nuclei and single hypernuclei. Furthermore, the effective field theory is designed 
for calculations at low and moderate nuclear densities. It is unclear how far the truncated 
theory can be extrapolated to high densities. We will leave these problems aside in the 
following and focus on the new effects arising from the flavor mixing. 
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Nuclear matter is bound over a wide range of X3 and xs- This can be deduced from 
Fig. |l|, where the binding energy is shown as a function of the baryon density. In part (a) we 
consider curves for various strangeness ratios at fixed isospin. At low and moderate densities 
an increase of the strangeness leads to higher energies. At high densities one observes the 
opposite trend; the curves become significantly softer for higher values of xs- Part {b) of 
Fig. indicates binding energies for various isospin ratios at fixed xs- Increasing the overall 
isospin always adds more repulsion to the system. The origin of this behavior is mainly a 
Fermi gas effect. The main contribution to the energy in Fig. |l| comes from the nonstrange 
baryons. Increasing the strangeness increases the fraction of the strange baryons leading to 
more symmetric systems with smaller energies. If the absolute value of the isospin increases 
the fraction of certain flavors increases leading to more asymmetric systems with higher 
energies. 

Typical flavor fractions are depicted in Fig. ^ In addition to the nucleons the A and E° 
fractions are nonzero at all densities for > and 2:3 7^ 0. The onset of the other flavors 
is determined by the condition 

vf> M*p . (74) 

The properties of the A and TP is governed by the flavor mixing. The effect is driven by 
the mean field of the p meson and does not occur in isospin saturated systems (0:3 = 0). As 
the lightest hyperon the A is expected to give the main contribution of the strange baryons 
at low densities. Remarkably, the baryon density of the TP is nonzero even if the condition 
Eq. (^) is not fulfilled. This is because the nuclear matter ground state is not a state filled 
with distinct A and S° flavors rather than a state filled with the mass eigenstates 1 and 2 
of the Dirac equations Eqs. ( |32D and (p^). 

The situation is illustrated in Fig. |^. The lower half of part (a) indicates the contributions 
of the mass eigenstates 1 and 2 to the thermodynamic potential according to Eq. The 
labels are chosen such that solution 1 reduces to the A and solution 2 to the E° in the limit 
Kn = 0. Below the density p\ only solution 1 gives a contribution. The mass eigenstate 
1 consists mainly of the A flavor with a small admixture of the T^ leading to nonzero S° 
density fraction in this region. It increases more rapidly when the second mass eigenstate 
emerges at p^. The system in part (b) of Fig. |^ has the same strangeness ratio as in part (a) 
but X3 = 0. As a consequence the mixing vanishes (V^ = 0). The lower half now indicates 
the thermodynamic potentials of pure A and T'^ flavors. The onset of the T'^ is governed 
by Eq. (|7^ which is not fulfilled for pb < p\' ■ The density and the contribution to the 
thermodynamic potential of the S° vanish in this region. 

The onset of the pure strange flavors depends strongly on the overall isospin and 
strangeness content of the system as indicated in Fig. ^ for the and Decreasing 
X3 delays the onset of the S+ and For sufficiently small values of X3 these two flavors 
disappear altogether. 

Another signature of A — flavor mixing is a nonvanishing value of the mixing density 
Eq. (|5lD . At high densities the mixing density is comparable to the baryon densities of the 
strange flavors as indicated in Fig. |^. The effect is stronger if either the overall isospin or the 
strangeness is increased. This can be studied in Fig. ^, which shows the mixing density for 
various isospin and strangeness ratios. The kinks at low densities are caused by the onset 
of the S~. At low and moderate densities, the mixing density is very small. 
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Generally, we observe a strong parameter dependence. First, the coupling (7^^ ^ind 
therefore the flavor mixing is smaller for the set Gl. Second, in the high density regime the 
set Gl and G2 predict significant different flavor fractions. This is a common observation in 
relativistic mean field models. Different models and parametrizations which are equivalent 
at low and moderate densities can produce significantly different results when extrapolated 
to regimes of high densities p8|. 



VI. NEUTRON STAR MATTER 

In this section we apply our model to describe dense matter in neutron stars. Recently, 
the study of strangeness in dense matter has received considerable attention due to the 
possibility of kaon condensation |]13[. However, the question of kaon condensation is quite 



delicate and very sensitive to the employed parameters and models [|T^,|l^. It is therefore 
important to study all the relevant many-body effects which have impact on the equation 
of state. 

Throughout the last section we assumed that the time scales are sufficiently short such 
that weak processes can be neglected. As a model for cold neutron star matter we consider 
baryons, electrons and muons in chemical equilibrium with respect to weak decays p4 



This is equivalent to introducing two chemical potentials characterizing the conservation of 
baryon number and electric charge. This situation is realized by adding contributions of free, 
relativistic electrons and muons to the baryonic equation of state. The chemical potential 
for the electric charge is introduced by adding 

= Mq - ^e7°^e - ^^7°^m) , (75) 

to the effective Lagrangian in Eq. (]^). The charge density of the baryons is given by 

J°Q = Tr(fi7°[Q,5]) , (76) 
where the charge matrix is given in terms the isospin and hypercharge operator by 

Q = \y + T,. (77) 
Overall charge neutrality is achieved by imposing 

< Gl J°Q - ^e7°^e - ^^.T^^^IG >= , (78) 

on the ground-state expectation value of the charge density. The resulting equation of state 
for the two parameter sets Gl and G2 are shown in Fig. |^. At low and moderate densities 
the two parameter sets are essentially equivalent. At higher densities the set G2 produces a 
significantly softer equation of state. 

The corresponding density fractions are indicated in Fig. |^. At low densities the proton 
fraction is very small and the system resembles pure neutron matter. At Pb/p% ~ 0.75 the 
muons start to emerge. The first strange baryon, the S~, starts at Pb/p% ~ 2 followed 
closely by the A. Due to the flavor mixing there is also a small fraction of the S". The onset 
of the S~ is at ps/p^ ^ 4.5. As the last flavor, the S"*" emerges at Pb/p% ~ 7.5. The fraction 
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of the S° becomes noticeable around Pb/Pb ~ 5.5 where the second mass eigenstate starts 
to contribute to the equation of state as indicated by the black dot. In the usual scenario 
T6| , |4^ without the flavor mixing the S° would not be present in the system below this point. 



Also included in Fig. y is the mixing density. At high densities its size is comparable to the 
other strange flavors indicating strong flavor mixing. A more quantitative estimate of the 
flavor mixing at low densities can be obtained from Fig. Ij. Similar as in the previous figure, 
the onset of the second mass eigenstate is indicated by the black dots. Below that point the 
sum of the flavor densities equals the density of state 1, according to Eq. ([56|). It mainly 
consists of the A flavor, the admixture of the E° is very small ( < 0.1%). 



VII. FLAVOR OSCILLATIONS 



As discussed in section |III B| the A and S*^ flavor eigenstates are not the mass or energy 



eigenstates of the Lagrangian in Eq. ( PU[ ) rather than a superposition of the actual mass 
eigenstates. The true, time independent nuclear matter ground state arises as a state filled 
with particles corresponding to the mass eigenstates. Adding a distinct flavor to the ground 
state creates a perturbation which evolves nontrivially in time leading to flavor oscillations. 
To study this phenomenon in more detail we will restrict our considerations to the limit 
= = M*. In this case the algebraic structure simplifies considerably allowing us to 
give explicit analytic expressions yet preserving the main physical content. 

According to section [111 B| the general solution of the Dirac equations Eqs. (|32| ) and ( |33D 



which uncouples the Lagrangian Eq. (pOf) can be written as 



^'a = cos(a)^i + sin(a)^2 , (79) 
^'so = - sin(a)^i + cos(a)^'2 • (80) 



The fields \Efi and \E'2 describe the normal modes which can be quantized canonically 

= e-^^» E (^:(p)«(^nP)e-^^^ + rf:t(p)^W(p)e^P-'^ , (81) 



with E = + M*2 and with the potentials introduced in Eqs. ( P8| ) and (^). The 

operators h\ and d\ satisfy the usual anti-commutation relations. In the mean field approx- 
imation the nuclear matter ground state \G > contains positive-energy levels of particles 1 
and 2 filled to the chemical potentials (Fermi energies) 

Pi = cos^{a)p\ + sin^(a;)/iso , (82) 
P2 = sin^(a)/iA + cos^(a)/xso . (83) 

A state of flavor A with momentum p and spin s relative to the ground state can be con- 
structed by adding a particle 1 and 2 above the Fermi energies 

|y9A(p, s) >= (^cos{a)bl\p) + sin(a)6f (p)) \G > with E + V^ > pi . (84) 

The normalization is chosen such that at = 
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J^d'x < MP,s)\^Al'^A\Mks') > ={27rf—6ss'6^'\p-k){l + N^) , (85) 
J^d'x < Mp, s)l^so7°^so|</^A(i, s') > = (27r)3 - k)N^° , (86) 



with 

J^d^x < G|^a7°^a|G' >= and J^d^x < G|^207°^20 >= A^^" , 

for the number of A and S*^ in the ground state at xq = 0. At a later time xq > the 
transition amphtude of finding the system in the A flavor state is given by 



dpo 



(87) 



(2n? 6ss'S'^'\p - k) /^e-^^'»"«u('^)(p)7°G^^(p)7%(^)(p) 



The transition amplitude involves the AA-Greens function which can be computed by using 
the representation Eqs. (^) and (^ for the field operators ||4^ 



G^\x -y) = < G\T[m^{x)m^{y)]\G> _ (88) 

= cos2(a) < G|r[^i(x)^i(i/)]|G > +sin2(a) < G\T[m ^{x)^ 2{y)]\G > 
= cos^{a)G^{x — y) + sin^(a;)G'^(s — y) . 

The momentum-space representation of the Greens function of particles 1 and 2 is 



G\p)=t 



li^- V.+ M*)5ipo - - E)e(/x. - V,' - E) . (89) 



{p - 1^)2 - M*2 + ie E 



After performing the po integration in Eq. { ^7\ ) the transition probability of finding the 

system in the state \<fA{p, s) > at time Xq follows to 



E 



(90) 



M 



— (^1 — sin^(2a;) sin^ 



(V° - V,') 



-Xo 



BiE + Vl> - /ii)e(E + - /is) . 



Similar as in Eq. (|8^ one can also define a S*^ state 

|v3so(p,s) >= (sin(a)6^"^(p) - cos(a)6f (p)) |G 



> 



(91) 



The probability for a transition of the state \ip\ > to the state ly^^o > involves the mixed 
AE°-Greens function. Repeating the steps which lead to Eq. (RO), it becomes 



E' 



*2 



sin^(2Q;) sin^ 



-Xq 



e{E + V^ - fXi)e{E + V^ - fX2) ■ (92) 



^We consider the system in a finite volume V and take the thermodynamic limit at the end of the 
calculation. 
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Thus the probabihties for finding the system in the A or S° fiavor state oscillate with 
frequency 

= ^ , (93) 

arising from the interference of particle 1 and 2. The nuclear matter ground state also offers 
the possibility to create A hole states by destroying a particle 1 and 2 below the Fermi 
energies 



liflip, s) >= [cos{a)bl{p) + sm{a)bl{p)) \G > with E + V^ K^a . (94) 
The probability of finding the system in the A hole state at a later time Xq 

P,Tixo) = - sin^(2«) sin^[ ^^' -/'K o])e{f^, -E- V^)Q{^,, -E- V^) , 

is identical to the expression for P^'^ in Eq. (|90[) except for the arguments of the step 
functions. 

The fiavor oscillations are governed by the frequencies ujp and by the factor sin(2a). The 
system unmixes if sin(2Q;) = 0, maximal mixing occurs for sin(2Q;) = 1 which corresponds to 
sin(a) = cos(a) = \/2/2. For neutron star matter these quantities are indicated in Fig. |^. 
The results were obtained by using M* = (M^ + M^)/2 0. The frequencies in the upper 
part of Fig. ^ steadily increase with increasing density. The factor sin (2a) is shown in the 
lower part. The mixing angle in Eq. (^) is determined by the mean field of the p meson 
and by the difference — V^o ■ The mixing is large at low densities because here the matter 
is very neutron rich. The system becomes more symmetric with increasing density leading 
to smaller mixing angles. At the minimum {pb/p% ~ 2) the S~ starts to emerge as the first 
strange baryon and the mixing angle increases again. According to our remarks at the end 
of Sec. we observe a sizable parameter dependence. The set G2 predicts higher mixing 
frequencies and stronger mixing. To obtain an estimate for the size of the effect, the results 
in the lower part of Fig. ^ can be compared with the rather small value sin(2Q;"'^) f» 0.02 
characterizing the mixing in the vacuum induced by the mass difference of the up and down 
quarks [|^. Thus, in the medium flavor mixing is much larger than typical isospin violating 
effects and thus might be easier to observe. 



VIII. SUMMARY 

In this paper we study strange nuclear matter based on concepts of effective fleld the- 
ory. Our starting point is a chiral effective Lagrangian containing the baryon octet, the 
Goldstone boson octet, the vector meson octet and a light scalar singlet. Guiding principle 
are the underlying symmetries of QCD which, in principle, implies that an inflnite number 
of interaction terms has to be considered. We truncate the effective theory based on the 



Numerically, this is in fact a good approximation for calculating the equation of state. 
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observation that at low and moderate densities the meson mean fields are small compared to 
the baryon masses and thus provide useful expansion parameters. We develop a mean field 
model based on the one-loop approximation of the thermodynamic potential. It contains 
the one-body contributions of the baryons and a nonlinear potential which is a fourth order 
polynomial in the meson mean fields. These fields are interpreted as relativistic Kohn-Sham 
potentials and nonlinear interactions between the meson fields parametrize density depen- 
dence which, in principle, is beyond the one-loop or mean field level. A necessary condition 
is that the model parameters can be accurately calibrated to observed properties of ordinary 
nuclei and hypernuclei. 

In the nonstrange sector we use the coupling constants which have been obtained in an 
EFT description for normal nuclei fl^. The scalar hyperon couplings are determined by 



using phenomenological information on the hyperon potential in nuclear matter. However, 
the nonlinear meson-meson couplings involving the meson remain largely unconstrained. 

In spite of this difficulty the EFT description leads to new and interesting effects. The 
most important observation is that the D-type coupling between baryons and vector mesons 
leads to a nondiagonal vector self-energy in the A — S° sector of fiavor space. As a conse- 
quence A — S° fiavor mixing arises. Our basic goal is to provide a first orientation of the 
formalism by studying nuclear matter. Although, the discussion in nuclear matter is an 
oversimplification we believe it is useful for providing a concrete description and for exam- 
ining qualitative features of the fiavor mixing. We discuss various nuclear matter systems 
with different isospin and strangeness content as a simple model for multi-strange systems. 
We find that nuclear matter is in a state of mixed fiavor rather than in a state with dis- 
tinct A and S° particles. This implies that systems which contain A hyperons always have 
a small admixture of S° hyperons. At low and moderate densities the features are small 
but we expect that signatures will survive in calculations of very heavy and asymmetric 
hypernuclei. 

The model was then used to study dense matter in neutron stars. In the presence of 
strangeness the equation of state suffers considerably from model and parameter dependence. 
Particularly, predictions for the onset of kaon condensation depend sensitively on model 
features relating to the nucleon-hyperon and hyperon-hyperon interaction. It is therefore 
important to estimate the influence of all the relevant many-body effects on the high-density 
equation of state. We flnd qualitative new features for the occurrence of the individual flavors 
in a neutron star. Due to the flavor mixing the A and start to emerge at the same value 
of the baryon density in contrast to the usual scenarios were the onset of the S° is delayed. 

Similar to the phenomenon of neutrino oscillations, the time independent ground state 
of nuclear matter responds to external perturbations with flavor oscillations characterized 
by distinct frequencies. In contrast to the neutrino oscillations which primarily occur in the 
vacuum the A — S° mixing is a true many-body effect and vanishes at zero baryon density. 
Moreover, mixing of the A and E° can also arise in the vacuum if isospin violating effects are 
taken into account, however, we find that the effect in the medium is considerably larger. 

To summarize, concepts and methods of effective field theory which have been applied 
to ordinary nuclei can be extended to incorporate strangeness in nuclear structure. The 
complexity of the underlying SU{3) symmetry leads to a large number of interaction terms 
which are not well constrained in the strange sector. To achieve more predictive power 
future studies are needed to determine the coupling constants by incorporating the presently 
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available information on single and double hypernuclei. 
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TABLES 



TABLE L Chemical potentials {up = fip — ^) and effective baryon masses. 



F 








P 

n 

A 
SO 
S+ 

s- 


1 1 ^ _l_ 1 1 * 

+ 2/^3 

/V-B - 
PtB + IJ'S 
fJ'B+IJ-S 
I^B + IJ-3 + I^S 
I^B- IJ-3 + PS 
PB + + 
MB - 5/^3 + 2/i5 


„aj , ,0 _i_ „0 J, 1 1 „P „ 

9n^ +9n'Po+29nPo 
gffUj^ + g%(t)o - y^Po 

g^u^ + gUo 
ffs^° + 5s'^o + gj:Po 

9%^° + 9t(t>o - 9sPo 
g^Lo'> + gt<Po + 

g%,u^ + gUo-hy^ 


Mjv — 5^9? 

Mn - 9%^ 
Ma - 9i^ 
Mn - ghv 
Mn - gl:<p 
Mn - gh<f 
Mn - giv 
Mn - gsV 



TABLE IL Baryon and vector meson masses (in MeV). 



Baryons 


Vector mesons 


Mn = 939 


= 782 


Ma = 1115.6 


rUp = 770 


Ms = 1193 


= 1019 


Ms = 1315 





TABLE III. SU{3) relations for the relevant uj,p and (j) baryon coupling constants. 



Vertex 


Coupling constant 


NNlo 


5jv = gs cos{9) + ^ 


/lgFsm{0) - -^go sin(0) 


NNcj) 


gtr = 9S sm{e) - y 


l^gF cos{e) + ^go cos{9) 


NNp 


gN = V^igp + go) 


AAlo 


g% = gs cos{e) - y 


'Igo sva.{e) 


AAcj) 


gt = ss sm{e) + ^ 


Igo cos{9) 


SSo; 


g% = gs cos{9) + y 


'Igo sm{9) 


SS(/. 


gf. = gs sm{e) - ^J 


^go cos{9) 


SS/) 


5e = V^gp 




AS/9 


5aE = V 3f^ 






gs=gscos{e) - y 


'Igp sin(0) - -^go sm.{9) 




gt = gs sm{e) + ^ 


IgF cos{9) + -^go cos{e) 




5s = V2{gF - go) 
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TABLE IV. Parameter Sets. Also included is the value of the baryon density p% at nuclear 
matter equilibrium. 





Gl 


G2 




0.53963 


0.55410 


N / 


0.78532 


0.83522 




0.96512 


1.01560 


0.69844 


0.75467 




6.3415 


10.462 




-286.15 


21.359 




0.48322 


4.7310 




-64.952 


8.3851 




518.48 


430.26 




-1.8063 


2.7532 


p% [fm-^] 


0.153 


0.154 



FIGURES 

FIG. 1. Binding energy of nuclear matter, (a) at fixed isospin ratio for various strangeness 
ratios. The strangeness ratios are xs = 0.2,0.3,0.4,0.5 from the bottom to the top. (b) at fixed 
strangeness ratio for various isospin ratios. The isospin ratios are = —0.1, —0.15, —0.2, —0.25 
from the bottom to the top. The parameters axe based on set Gl. 

FIG. 2. Flavor fractions as a function of the total baryon density. The parameters are based 
on set G2. 

FIG. 3. (a) Flavor fractions and contributions of mass eigenstates 1 and 2 to the thermodynamic 
potential, (b) Flavor fractions and contributions of pure A and flavors to the thermodynamic 
potential. The system in part (b) exhibits no flavor mixing (x^ = Vm = 0). 

FIG. 4. Flavor fractions of the pure strange flavors at fixed strangeness ratio for various isospin 
ratios. The parameters are based on set Gl. 

FIG. 5. Mixing density for various strangeness and isospin ratios. The parameters are based 
on set G2. 

FIG. 6. Equation of state of neutron star matter for the two parameter sets Gl and G2. 

FIG. 7. Density fractions in neutron star matter for the various baryon and lepton flavors. The 
parameters are based on set G2. 
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FIG. 8. Density ratios of the mixed A and S° flavors in neutron star matter. 



FIG. 9. Flavor mixing frequencies and the factor sin(2a) in neutron star matter. 
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(b) 



FIGURE g 
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